Entangling power of quantum chaotic evolutions via operator entanglement 
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We study operator entanglement of the quantum chaotic evolutions. This study shows that prop- 
erties of the operator entanglement production are qualitatively similar to the properties reported 
in literature about the pure state entanglement production. This similarity establishes that the 
operator entanglement quantifies intrinsic entangling power of an operator. The term 'intrinsic' 
suggests that this measure is independent of any specific choice of initial states. 
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Entanglement is a unique quantum phenomenon which 
continues to baffle and surprise. Recently, this remark- 
able property of quantum mechanical systems has been 
identified as a quantum resource whose production is an 
elementary prerequisite for any quantum informational 
and computational tasks Q. This basic task is accom- 
plished by unitary transformations, i.e., a given unitary 
operator acts on a product state and transforms that 
state into an entangled state. Then this given operator 
is referred to as an entangling operator. However, entan- 
gling power of different unitary operators are naturally 
not the same. Two different unitary operators operat- 
ing on two identical initial states can produce different 
entanglement. 

Investigations of different quantum signatures, like 
spectral statistics of the quantum chaotic Hamiltonian, 
phase space scarring, fidelity decay, etc., of classically 
chaotic systems is the subject of "quantum chaos" |2|. 
Recent studies have shown that entanglement in chaotic 
systems can also be a good indicator of the regular to 
chaotic transition in its classical counterpart 0, 0, 0, 
ISHHSESIllElElEi- Some of these studies 
have reported that the presence of chaos enhances the 
entanglement production rate However, in our ear- 

lier works, we have observed saturation of entanglement 
production for strongly coupled strongly chaotic systems 
0, 0, 0- This saturation of the entanglement produc- 
tion is a statistical property which we have modeled by 
random matrix theory (RMT) [6| . Coupling strength be- 
tween two chaotic subsystems is another important pa- 
rameter for the entanglement production. For example, 
in case of weak coupling, the entanglement production 
is higher for sufficiently long time corresponding to non- 
chaotic cases H- 

Recently, rather than focusing on the entanglement 
evolution of single initial product states, the authors 
studied glob al entangling properties of coupled chaotic 
systems jl5| using the coupled kicked tops model 0, @ . 
Following Ref. [16j, they have considered different en- 
sembles of initial product states and have studied ensem- 
ble average of the entanglement production as a measure 
of entangling power of coupled kicked tops time evolu- 



tion operator Ut- Their results have shown that the 
entangling power of Ut is quantitatively, as well as qual- 
itatively different for two different ensembles of initial 
states, one that is averaged over all product coherent 
states, the other including all appropriate product states 
jl, r | . Thus this averaging gives a more global measure of 
entanglement. 

In the present paper, we study operator entanglement 
of quantum chaotic evolution operator. We find that the 
operator entanglement production behaves qualitatively 
similar to pure state entanglement production [8|- This 
similarity justifies that the operator entanglement is an 
intrinsic measure of the entangling power of a given oper- 
ator. The term 'intrinsic' is referred to the fact that this 
measure is independent of any specific choice of single ini- 
tial state or some ensemble of initial states. Indeed the 
operator entanglement is only a property of the operator 
and makes no reference to states it acts on. Further- 
more, it has been shown in Ref. that the operator 
entanglement is related to the measure proposed in |l6| . 
However, the operator entanglement is easier to calculate 
than the measure proposed in |16| . 

For a formal definition of this measure, let us start with 
a simple algebraic transformation on matrices called ma- 
trix reshaping |l8j . Consider a rectangular matrix A with 
elements Ay , i = 1, . . . , K and j — 1, . . . , L. This ma- 
trix can be reshaped into an one-dimensional vector \A) 
by putting its elements row after row into lexicographical 
order of size KL, i.e. 



' m\A ) — Aij where m = (i — 1)L + j, 
i=l,...,K; j = 1, . . . ,L. 



(1) 



Following is a very simple example of the matrix reshap- 
ing [3 : 



.4 



An A 12 
A21 A22 



<-» \A) ={4u A 12 A 31 A 22 } T . (2) 



Vector \A) can be considered as an element of operator 
Hilbcrt space or Hilbert-Schmidt (HS) space. The scalar 
product between any two elements (A, B) of a HS space 
Hhs is defined as (A\B) = TiA^B. Consequently, the 
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HS norm of a matrix is equal to just the norm of the 
associated vector, i.e. ||A||!f g = = J2 m \ a m\ 2 - 

Let us now consider an arbitrary unitary operator U 
operating on a bipartite state space Ji = Hi ®Ti 2 , where 
dimWi = N < dimH 2 = M, and U € We 
are interested to measure entangling power of the uni- 
tary operator U. This unitary operator can be expanded 
in terms of complete orthonormal operator basis states 
{\A m ) (g> \B a )} as 



\U] 



N 2 M 



y j y 1 x ma 



\A 



\B n 



(3) 



where \U), \A m ) and \B a ) are the associated reshaped 
vectors of the matrices U, A m and B a . The operator ba- 
sis states ||A m )} and {|£? Q )} are orthonormal in the 
sense that they satisfy (A m \A n ) = Tr (A^An) = 5 mn 
and ( B a \Bp ) = Tr(i?£ Bp) = S a/3 . A very simple 
example of a complete orthonormal operator basis is 
{J2, <Ji}/y/2, where I 2 is the 2x2 unit matrix and oVs 
are the Pauli spin matrices. Now \U) can be considered 
as a vector in the composite HS space, Hhs ® Whs- We 
can apply Schmidt decomposition to | U ) and get 



\U) = 



N 



Am, \A Tl 



\B„ 



(4) 



where {A m } are the singular values of the rectangular 
matrix X, and {|A m )} and {|_B m )| are the new or- 
thonormal basis states. {A m } can also be identified as 
the nonzero eigenvalues of operator reduced density ma- 
trices. The operator Schmidt decomposition has also 
been proposed by Nielsen et al |l9j. However, applica- 
tion of the matrix reshaping operation has put operator 
Schmidt decomposition and state Schmidt decomposition 

on a same footing. Here we notice (U\U) = J2m=i ^ m = 
Tr(WU) = NM, i.e. the vector \U) is not normalized. 
To normalize \U), we define A m = X m /NM. We now 
define von Neumann entropy Sy(U) and linear entropy 
Sl(U) of the operator entanglement respectively as 

N 2 N 2 . . 

S V (U) = - y X m In A m - - Y -22- In (5) 
v\ ) ^ m m MM NM v ; 



m— 1 

N 2 



s L (u) = i-X)^ = 1 -E 



\2 



N 2 M 2 ' 



(6) 



These measures have already been utilized to study en- 
tanglement capability of qudit gates ^3 • Both the mea- 
sures give qualitatively similar results, but the von Neu- 
mann entropy is a more acceptable measure. Therefore, 
in the present paper we prefer Sy(U) to investigate the 
entangling power of quantum chaotic unitary operators. 

We use coupled kicked tops as our model of coupled 
chaotic system 0,0- The time evolution operator, de- 
fined in between two consecutive kicks, corresponding to 



the coupled kicked tops is given by 

u T = m 2 (Ui ® u 2 ) = u{ 2 [(u?u() ® (ugui)] 

where the different terms are given by, 



(7) 



U 



exp y-i^Jvi) . u - 
U{ 2 = exp 



k 



exp 



2j t z ' 



(8) 



and i = 1,2 represents two different tops. The term Xj( 
describes free precession of the top around y axis with 
angular frequency tt/2, represents a torsion about z 
axis by an angle proportional to J z with the proportional- 
ity factor k/2j, and U{ 2 is the coupling between the tops 
using spin-spin interaction term with a coupling strength 
of e/V7tf2- 

Here we study the evolution of the operator entan- 
glement of Ut- More precisely, we study the operator 
entanglement of as a function of the time step n. In 
our calculation, we choose complete orthonormal opera- 
tor bases corresponding to each subsystem as, 



A n 



|mi)(ni| and Bp = \m 2 )(n 2 \ 



(9) 



where 



a = N(mi + ji) + (ni + ji + 1); (mi,ni) = -jx, 

(3 = M(m 2 +j 2 ) + (n 2 +j 2 + 1); (m 2 ,n 2 ) = -32, ■ ■ ■ ,h 

and N = 2ji + 1, M = 2j 2 + 1. Applying matrix reshap- 
ing operation to and to the orthonormal basis states, 
we can expand \Uj, ) as, 



N M 2 - 



U t) = 2^2^ u a0 {n)\A a )®\Bp). 



(10) 



= 1 /3=1 



Following the procedure discussed above, we investigate 
the operator von Neumann entropy Sy(U^) as a func- 
tion of n. This study will show how entangling power of 
chaotic evolution changes with time steps. 

In Fig^ we have presented our results for the opera- 
tor entanglement production in coupled kicked tops for 
the spin j\ = j 2 = j (say) =10. As we go from top 
to bottom windows, coupling strength is increasing by 
a factor of ten (e = 10~ 4 to e = 1.0). For each cou- 
pling strength, we have studied operator entanglement 
production for four different single top kicked strengths 
k, whose corresponding classical phase space picture has 
been presented in Rcf. 8] . Here we are presenting a brief 
qualitative description of those phase space pictures to 
correlate the effect of underlying classical dynamics on 
the operator entanglement production. For k = 1.0, the 
phase space was mostly covered by regular orbits, with- 
out any visible stochastic region. For k = 2.0, most of the 
phase space was covered by the regular region, but with 



3 




O 200 400 600 800 1000 




O 200 400 600 800 1000 



I 




O 20 40 60 80 100 



FIG. 1: Time evolution of the operator von Neumann en- 
tropy corresponding to the coupled kicked tops is presented 
for different coupling strengths and for different underlying 
classical dynamics. Solid line represents k = 1.0, dotted line 
corresponds to k = 2.0, dashed line is for k — 3.0 and dash-dot 
line represents k — 6.0. 

a thin stochastic layer at the separatrix. For the change 
in the parameter value from k — 2.0 to k = 3.0, there 
was significant change in the phase space. At k = 3.0, 
the phase space was of truly mixed type. However, the 
size of the chaotic region was much larger than the reg- 
ular region. Finally, when k — 6.0, the phase space was 
mostly covered by the chaotic region, with few very tiny 
regular islands. 

Let us first discuss the case of weaker coupling e = 
10~ 3 , whose results are presented in the top most window 
of Fig^ Here we observe larger operator entanglement 
production for the nonchaotic cases than the chaotic one. 
Therefore, for very weak coupling, the presence of chaos 
actually suppresses the entangling power of U^. This 
property has already been observed in 0, El . Present 
observation shows that the suppression of entanglement 
by chaos is an intrinsic property of coupled chaotic sys- 
tems and is independent of any specific choice of initial 
state. 

The results corresponding to e = 10~ 2 are presented in 
the second window from the top of Fig^ In this case, at 
least for initial time steps (n < 100), we have observed 
suppression of operator entanglement by chaos. However, 
for larger n, the operator entanglement corresponding to 
both the nonchaotic cases (k = 1.0, 2.0) eventually satu- 
rates at a value lower than the values corresponding to 
the mixed (k — 3.0) and the chaotic (k = 6.0) cases. 
Here we also notice that the operator entanglement pro- 
duction corresponding to the chaotic case (k = 6.0) is 
always larger than the mixed case (k = 3.0) 

We now come to a reasonably stronger coupling 
strength e = 10 _1 , for which results are presented in 
the third window from the top of Fig^ Here again we 



observe saturation of the operator entanglement produc- 
tion for all the cases. However, the operator entangle- 
ment production rate is now much higher than the pre- 
vious cases. These saturation values are clearly differ- 
ent for nonchaotic, mixed and chaotic cases. For the 
chaotic case k — 6.0, our numerical estimation shows 
that the saturation value is approximately ln(0.6./V 2 ) = 
ln(0.6 x 441) ~ 5.57. On the other hand, the satura- 
tion value corresponding to the mixed case k — 3.0 is 
slightly lower than this value. But the saturation val- 
ues corresponding to the nonchaotic cases are distinctly 
lower than the other two cases. 

Finally, we discuss the case of very strong coupling 
(e = 1.0) for which results are presented in the bottom 
window of Figd Here, due to the strong coupling, the 
over all coupled system is chaotic irrespective of underly- 
ing classical dynamics of the individual subsystems. Con- 
sequently the saturation values of the operator entangle- 
ment production are almost equal to ln(0.6iV 2 ) ~ 5.57 
for all different values of k. 

An important outcome of this study is the observation 
of operator entanglement saturation at around ln(0.67V 2 ) 
in case of chaotic subsystems which are coupled very 
strongly. Consequently, from our previous knowledge of 
the saturation of the pure state entanglement , we ex- 
pect that the distribution of the eigenvalues of the oper- 
ator reduced density matrices (RDMs) will follow RMT 
prediction. In Fig[2 we have presented the distribution 
of the eigenvalues of the operator RDM corresponding to 
coupled kicked tops for different Hilbert space dimension 
M of the second top. Here we have fixed the dimension 
of the Hilbert space of the first top at N = 2j\ + 1 = 21. 
The solid curve is representing the RMT predicted La- 
guerre distribution, i.e., 

, m _ N 2 Q V(A ma x-A)(A-A mm ) 

2tt A ^ 

xmax _ J_ ( i , J_ _i_ 2 \ 
min N 2 \ Q ^Qj 

where A € [A mm , A ma x] and Q = M 2 /N 2 ; and the his- 
tograms are the numerical results obtain from coupled 
kicked tops model. This figure shows clear agreement 
between RMT prediction and numerical results. Using 
above distribution and following the procedure presented 
in 0, we can analytically estimate the operator entan- 
glement saturation for different dimension ratios Q. This 
calculation is straightforward, therefore we do not pursue 
this further here. 

Above studies on the operator entanglement produc- 
tion show some qualitative properties which are common 
to the observed properties of the pure state entangle- 
ment production (see 001 El)- For instance, (1) the 
operator entanglement, in general, supports 'more chaos 
more entanglement hypothesis; (2) the operator entan- 
glement production has also one statistical upper bound 
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FIG. 2: Distribution of the eigenvalues of the operator RDMs 
of coupled kicked tops (k — 6.0, e — 1.0). Time step n has 
been taken in the saturation region. TV = 2ji + 1 = 21. Solid 
curves correspond to the theoretical RMT predicted distribu- 
tion function. 

which can be explained by RMT, and (3) for the weakly 
coupled cases, chaos suppresses the operator entangle- 
ment production. These similarities confirm that these 
are generic intrinsic properties of the entanglement pro- 
duction under chaotic evolution. 

Eigenstates of any entangling operator also show some 
intrinsic entanglement properties of that operator. How- 
ever, in many cases, eigenstates fail to give any clue of 
the entangling power of the operator. Consider an en- 
tangling operator U = exp(— iaJ Zl <8> J Z2 ) where J Zi 's 
are usual angular momentum operators. Its entangling 
power is determined by the parameter a. For a definite 
a, this operator can create a maximally entangled state 
in case of a very special initial state. The eigenstates of 
J Zi 's, i.e. {|mi,m 2 )}, are also the eigenstates of U with 
eigenvalues {exp(— iamirri2)} ■ These eigenstates are all 
uncntangled, but many of them are degenerate. Any lin- 
ear superposition of degenerate eigenstates can form en- 
tangled eigenstates, which are independent of the param- 
eter a. For example, \(j>) = (\m 1: m 2 ) + 1 — m 1; — m 2 ))/V2 
are the entangled eigenstates of U for mi = or m 2 = 0, 
but TOi,m2 are not simultaneously equal to zero. These 
eigenstates belong to the largest degenerate subspace of 
dimension Aj + 1 with eigenvalues of unity. However, 
the maximum possible von Neumann entropy of these 
entangled eigenstates can be In 2. This is a maximally 
entangled state only for the spin j = 1/2. Besides, there 
are some degenerate subspaces with dimensions less than 
Aj + 1, in which we can construct eigenstates with von 
Neumann entropy larger than In 2 but much less than 
the maximum possible value ln(2j + 1). Above all en- 
tanglement of these eigenstates are independent of a. 
On the other hand, operator entanglement of U is a 
function of a and it shows qualitatively similar behav- 



ior to the pure state entanglement production by U. In 
case of the operator U^, we have seen that its entan- 
gling power increases with n. But the eigenstates of 
are independent of n and consequently it is impossible 
to distinguish the entangling power of Uj, for different 
n. Let us consider an interesting unitary operator U p — 
cxp(— ipJ Zl ) ® exp(— ipJ Z2 ) — exp[— ip(J Zl + J Z2 )] whose 
eigenstates are again { |mi , 7712) } . This operator is clearly 
not an entangling operator. However, due to degeneracy, 
there are some eigenstates like \<j) p ) = N~ 2 J2 m \m, — m) 
are maximally entangled. Moreover, from different de- 
generate subspaces, we can construct eigenstates with 
von Neumann entropies In A where X can be any real 
number between 2j + 1 to 1. Therefore, in this case, en- 
tanglement of the eigenstates gives an incorrect estima- 
tion of the entangling power of U p . However, operator 
entanglement of U p is equal to zero. 

In summary, we have studied operator entanglement 
production of the chaotic evolutions. This study shows 
that the behaviors of the operator entanglement produc- 
tion are qualitatively similar to the reported behaviors 
of the pure state entanglement production by the chaotic 
evolution operators. Consequently, this study establishes 
that the operator entanglement is an intrinsic measure of 
the entangling power of an operator. Moreover, we have 
pointed out that, when the eigenstates of a given operator 
fail to give any clue of its entangling power, the operator 
entanglement can give some estimation. As fur as RMT 
is concerned, we have demonstrated another realization 
of the Laguerre distribution. 



* Electronic address: jayendra@prl.ernet. in 
t Electronic address: arul@physics.iitm.ac.in 
[1] M. A. Nielsen and I. L.Chuang, Quantum Computing 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 
[2] F. Haake, Quantum Signatures of Chaos, 2nd ed. 

(Springer- Verlag, Berlin, 2000). 
[3] K. Furuya et ai, Phys. Rev. Lett. 80, 5524 (1998). 
[4] P. A. Miller and S. Sarkar, Phys. Rev. E 60, 1542 (1999). 
[5] A. Lakshminarayan, Phys. Rev. E 64, 036207 (2001). 
[6] J. N. Bandyopadhyay and A. Lakshminarayan, Phys. 

Rev. Lett. 89, 060402 (2002). 
[7] A. Tanaka et al, Phys. Rev. E 66, 045201 (R) (2002) ; H. 

Fujisaki et ai, Phys. Rev. E 67, 066201 (2003). 
[8] J. N. Bandyopadhyay and A. Lakshminarayan, Phys. 

Rev. E 69, 016201 (2004). 
[9] A. Lahiri and S. Nag, Phys. Lett. A 318, 6 (2003). 
[10] A. J. Scott and C. M. Caves, J. Phys. A 36, 9553 (2003). 
[11] Ph. Jacquod, Phys. Rev. Lett. 92, 150403 (2004). 
[12] X. Wang et al, Phys. Rev. E 70, 016217 (2004). 
[13] S. Chose and B. C. Sanders, Phys. Rev. A 70, 062315 
(2004). 

[14] J. N. Bandyopadhyay and A. Lakshminarayan, Pramana 

64, 577 (2005) ; e-print [quant-ph/050m3] 
[15] R. Demkowicz-Dobrzariski and M. Kus, Phys. Rev. E 70, 



5 



066216 (2004). 

[16] P. Zanardi et al, Phys. Rev. A 62, 030301(R) (2000). 

[17] X. Wang et al, Phys. Rev. A 67, 042323 (2003). 

[18] K. Zyczkowski and I. Bengtsson, Open Syst. Inf. Dyn. 



11, 3 (2004) ; e-print quant-ph/0401119 
[19] M. A. Nielsen et al, Phys. Rev. A 67, 052301 (2003). 



